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SUMMARY 



This is an investigation of the equations of motion and 
physical parameters involved in stabilizing the initial flight 
of a vertically launched rocket by means of a booster rocket 
pin-connected below the main rocket. The system is designed 
to stabilize the flight in its early stage before the aerodynamic 
control surfaces become effective. Stability of the system is 
dependent on the pendulum action of the booster rocket. 

The equations of motion were derived from Lagrange's 
generalized momentum equation. The differential equations 
thus obtained were not solved but were tested for stability by 
means of Routh's stability criteria. The ratio of the mass of 
the main rocket, M^, to the mass of the booster rocket, M^, 



Ml 



Ml 



Ml 



The system involving a mass ratio- — = 7.75 was found 

Mp 



to be unstable under all conditions. However, the system in- 
Mi 

volving a mass ratio = 1.5 was determined to be stable in 
the range 1.62 ■< 4.54 x 10^®, where ^ is defined as the ratio 

of the distance ^ from the center of gravity of the booster M 2 
to the pin connecting the strut to the main rocket M^, divided 
by the radius of gyration, k 2 i of the booster M^. In this range, 
for any given value of , stability was uniquely determined by 



one value of the ratio a = ^ , where /, is the length of the 
strut fronn the main rocket to the booster rocket ^ 2 * Thus, 
for a given booster, stability is primarily a function of the ratio 
^ , and for any given is uniquely determined. 

Although the system was found to be theoretically stable 

Ml ^ 

for the mass ratio--— = 1 .5 , the ratio turned out to be of 
M2 y*. 

such great magnitude as to make the system entirely impracti- 
cal for this particular mass ratio. 
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I. INTRODUCTION 

When a rocket or guided missile is initially launched there 
is a short period when its flight is unstable due to the fact that at 
low velocities the aerodynamic control surfaces are ineffective. 

Various procedures have been adopted to stabilize this 
Initial part of the flight. The most obvious and most successful 
of these has been the use of guide rails, or launching towers which 
hold the missile on its course until it has attained a velocity at 
which the control surfaces become effective. But this method 
requires a bulky launcher, and special care is necessary to main- 
tain the rails in proper alignment. 

A modification of the guide rail system is the "short- 
length" launcher which is really a guide rail whose length has 
been reduced to a minimum by greatly increasing the initial ac- 
celeration and thus the velocity of the missile. This increased 
acceleration presents a great many problems in component design 
to resist the tremendous accelerative forces. 

During World War II the Germans were fairly successful 
with a unique approach to the problem in their well-known "V-2" 
rocket. (Reference 1). The "V-2" was launched from a near- 
vertical position with no external restrictions. Instead it em- 
ployed four carbon vanes mounted in the jet stream and activated 
by the auto-pilot to maintain the missile on its proper flight path. 
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The great difficulty here was that the vanes steadily burned out 
with a consequent diminishing of control. 

The subject of this thesis is a third means of initial flight 
stabilization based on the concept of the action of a pendulum and 
independent of both external guides and internal control mechan- 
isms . 

If the missile booster is constructed as a separate unit 
and attached to the main rocket by a pin connection a certain dis- 
tance below its own center of gravity it can be seen from Figure 
1 that any motion of the main rocket will impart a proportionate 
motion to the booster. Specifically, if the main rocket rotates 
about its center of gravity in a clockwise direction, the connec- 
ting strut will cause the booster to rotate also about its own center 
of gravity in a clockwise direction. Thus, if the main rocket 
turns off course in a clockwise direction the booster immediately 
turns in a clockwise direction also, that is, the booster tends to 
re-align its thrust with the axis of the main rocket. If, by a 
proper choice of moment arms, the booster can be made to over- 
shoot this position of thrust alignment slightly the thrust will have 
a horizontal component which acts on the strut to the main rocket 
and causes the main rocket to turn about its center of gravity in 
a counter-clockwise direction. Thus, the booster responds to an 



error and its response tends to counteract the error. Such a sys- 



. i 
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texn should produce a sinusoidal flight path if the proper moment 
arms and damping forces were applied. 

In this analysis frictional and aerodynamic forces were ne- 
glected to simplify the equations. The only damping force consid- 
ered (other than the inertia forces of the system) was the jet damp- 
ing force or the resistance of a jet stream to rotation. The other 
external forces considered were the weights of the main rocket 
and the booster, and the thrust of the booster. (The motor of the 
main rocket is considered not to be operating in this analysis). 
Taking into account only these forces the writer derived the equa- 
tions of motion of the system in two dimensions by means of La- 
grange's generalized momentum equation. Because of the com- 
plicated nature of these equations no solution of them was attempted. 
Instead the conditions for stability were investigated by means of 
Routh's stability criteria. 

A hypothetical system consisting of a five-second booster 
pin-connected below a "V-2" rocket was first investigated. For 
this system the ratio of the mass of the main rocket Mj, to the 

Ml 

mass of the booster M,, is = 7,75. To reduce the work of 

computation the lateral motion of the center of gravity of the sys- 
tem was set equal to zero and only the vertical motion of the 
center of gravity and the rotation of the main rocket and the 



booster were considered. 
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Since no stable solution was found, a system with mass 
Ml 

ratio !! 1.5 was next investigated, again neglecting the lateral 

Mg 

motion of the center of gravity of the system. In this case a range 
of stability was determined but it was of such a sensitive nature 
as to make impractical further investigation for this mass ratio 
in the general case with lateral motion taken into account. 
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II. DERIVATION OF THE EQUATIONS OF MOTION ' 

In deriving the equations o£ motion, only four external 
forces are considered, namely, the weight, Mj^g, of the main 
rocket, the weight, M^g* of the booster, the thrust, F, of the 
booster, and the jet damping force, D, of the booster. These 
forces are indicated in Figure 1. 



Figure 1. External Forces on Main Rocket and Booster. 



D 




F 
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The jet damping force, D, is a force which resists the ro- 
tation of the booster about its own center of gravity. If the booster 
is of mass then the rate 

of flow through the nozzle may 
dMx 

be written as , or As 

the booster rotates about its cen- 
ter of gravity the tangential ve- 
locity at the nozzle exit will be 
or 1^ . Therefore 

the damping force, D, will be 

equal to the rate of change of momentum, or M 2 * This 

t- XV T . i ml 

expression has the dimensions “y* 'L* ~ , or -:jr" • 

Frictional forces and aerodynamic forces are neglected in 

this analysis. 

Both the main rocket, M^, and the booster, M 2 * as- 
signed three degrees of freedom, namely, lateral, vertical, and 
rotational motion. Note, also, that only the two dimensional case 
is being treated. The coordinates of masses Mj and M^ are indi- 
cated in Figure 2. The connecting strut between M^ and M 2 is 
sumed to be of zero mass and of infinite stiffhess for purposes of 
this analysis, and consequently the coordinates of the center of 
gravity of the system are indicated on the straight line joining 
Mj and M£ in Figure 2. 





J 
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Figure 2. Schematic Diagram of Coordinate System and 
Dimensions » 



It is assumed that the propulsion system of the main rocket, 
Mj, will not be operating during the initial stages of launching and 
consequently Mj is a constant for this analysis. The booster, 
will of course be burning with a consequent change of mass during 
launching, and the rate of change of mass M2, say M2, will be 
considered, as will be the change of position relative to Mj and M2 
of the center of gravity of the system. 

The first step of the analysis is to eliminate the extraneous 
coordinates q^ , q^, q^, and q^ in order to write the equations of 
motion in terms of the rectilinear coordinates of the center of grav- 
ity of the system, q^, and q2» and the angular coordinates q3 and q^ 
of Mj and about vertical axes through their respective centers 
of gravity. (Since only the onset of instability is of interest, the 
analysis will be restricted to consideration of small variations of 
the angles q3 and q^.) Once this step has been completed, La- 
grange's generalized momentum equation will be used to deter- 
mine the equations of motion. 

ELIMINATION OF EXTRANEOUS COORDINATES 

The coordinates qj , q^ , , and qg are always redundant 

because of the rigid links in the system and the known position of 
the center of gravity of the two masses. Therefore these extran- 
eous coordinates may be eliminated by considering the constraints 




I 
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elaborated below. 

The first relation between the coordinatv^s of the individual 
centers of gravity and the center of gravity of the system may be 
found by taking moments about the center of gravity of the system 
indicated in Figure 3. 




Figure 3. Relationship of Mj and to the center of gravity of 
the system. 



Thus 9 



or 









( 1 ) 



f 
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Furtherixiore the horizontal and vertical separations of 
the centers of gravity of Mj and M 2 in Figure 4 are fixed by the 
length of the line, ^ , joining the two, and by the angle 
which makes with the vertical. 




Figure 4 . Relationship of and M 2 to the connecting strut. 
In Figure 4 the vertical separation is 






or V = P A 



and the horizontal separation is 

J - A-fz.-. 



( 2 ) 



or 



( 3 ) 
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A fourth equation follows from a consideration of the 
similar triangles indicated in Figure 3, 

-p r -A 

and, since the denominators of both sides of this equation are 
equal from Equation (1), It may be written 

A/'/i-f,) - <■*' 

Equations (1), (2), (3), and (4) provide four linear rela- 
tions between the coordinates q • q • 

q^ , and q^ , Consequently the former may be expressed in terms 
of the latter and thereby eliminated from" subsequent equations. 
Solving these equations the following expressions are obtained: 

From Equation (2) fc ~ ^ ( 5 ) 

From Equation (3) pa ~ ^) -t 

or ^ (6) 



Substituting Equation (6) into Equation (1) 
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^ /%) ~ y^ip + /% 

^ = / 

Substituting Equation (5) into Equation (4) 
/%[^y (fz -^r) 

/> /^ Af,y/^z / ^ 

Substituting Equation (7) into Equation (6) 

fr ^ r * -yCf.-i) '-^ 

/? = / ' -') 

^ = / 

From Equation (2) 

^ ~ ^ /> 

Substituting Equation (8) into this equation 



( 7 ) 



(8) 



( 9 ) 



V 



/f " * Ti "-<7^ -“O' 

f, ^ f, -> f'~ 707iJ 

^ f- 70»: ‘-■‘Cp-n , 10 ) 

These expressions for the extraneous coordinates in Squa- 
tions (7). (8), (9)> &nd (10) involve, in addition to the principal co> 
ordinates, the variables the distance between the centers of 

gravity of the two components, and ^ , the angle between the strut 
, and the line ^ , joining the centers of gravity of the two 
components. These two variables will now also be expressed in 
terms of the principal coordinates. 

In Figure 4 by the law of cosines 

^ ( 11 ) 

and by the law of sines 
A ^ Av 



whe nee 



( 12 ) 
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These relations are unfortunately complicated. There- 
fore, to simplify them for later use, the fact that the variations 
in and are small will be taken into account. Assuming 
these angles sufficiently small such that terms of order q* 
may be neglected, the following approximations can be made 

From Equation (11) 






From Equation (12) 








and since 










Thus the appropriate approximations are 



-fj - ^ 




> ( 13 ) 
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Substitution of the results of Equation (13) into Equations 
(7), (8), (9)f and (10) yields the following approximate relations 
between the extraneous and the principal coordinates 






( 14 ) 






fz 






( 15 ) 









(16) 



f' ~ ( 17 ) 

To derive the equations of motion by the method of La- 
grange the expression for the kinetic energy, T, must first be 
found for use in the Lagrangian equation 



-W r j _ ^ r 

J 






(18) 



where is the "generalized force" and q^. is • (Reference 2) 



The generalized force is employed rather than the potential 
energy because this is a non-conservative system. 

If ly denotes the moment of inertia of the mass about 
its own center of gravity, the kinetic energy of translation and 
rotation for the system may be written 
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( 19 ) 



DERIVATION OF EQUATION OF MOTION 
FOR COORDINATE qj 



In applying Equation (18) to the direction of the coordinate 
q. , the value of -4^ is first computed using the kinetic energy, 

T , of Equation (19) . 

^ ■!■ /% fr 4^ 

The extraneous coordinates q^ and q^ may be eliminated from 
this expression by differentiating Equations (14) and (16) respec- 
tively. From Equation (14) 



^ =/ 



and from Equation (16) 



f? t' ~ ''C fJ 

iZL. = / 

Therefore 



/v; 



^ -A 



^Xr- 






Af* 












Af. 
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The first term in the Lagrangian expression is then 






M, 



f^lft f*)l 



( 20 ) 

By inspection of Equation (19) t the second term of the La- 
grangian expression is 

^ = <? (21) 

The generalized force is most easily computed by con- 
sidering the work done by the external forces as the coordinate q. 
is varied, all other coordinates remaining fixed. 

In general, Work =^. A fj , or 
In particular, the force becomes, referring to Figure 2, 

^ Work _ ^"7^*' -T- D ^y) A 



B 
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But within the present approxinnations 






and ^ nf / 



hence F -t £> 



( 22 ) 



Substituting the results of Equations (20), (21), and (22) 
into Equation (18), the equation of motion in the q direction be- 
comes 



DERIVATION OF ECU AT ION OF MOTION 
FOR COORDINATE ^ 2 . 

The equation of motion in the q^ direction is obtained in 
an identical manner . 

From Equation (19) 







(23) 






and, differentiating Equation (15) 




Consequently 






/ 



Similarly, differentiating Equation (17) 
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so that 

t-' 

Therefore 

and the first term of the Lagrangian expression becomes 






Furthermore ^ (25) 

The generalized force t referring to Figure 2, 

^ = Work . 



and since 






- 



-*■ ~ -^ fy 



(26) 



ji 



II 

1 



I 

I 



I 
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Substituting the results of Equations (24), (25), and (26) 
into Equation (18), the equation of inotion in the direction be- 
comes 






( 27 ) 



DERIVATION OF EQUATION OF MOTION 
FOR COORDINATE q^ 



The first term of the Lagrangian expression is obtained in 
an identical manner to that employed for coordinates and . 
From Equation (19) 




and, differentiating Equation (14) 






Similarly, differentiating Equation (16) 









S! 



i 

1 

I 

i 

I 



■ 

( 

I 

I 

I 

1 

I 
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jr =. X.f-, 

0^ 

or 

^ =^/. ' ?^^('<f’-^^r^) 

and the first term of the Lagrangian expression becomes 



Mif\ ■^'f-'& 

From Equation (19) 



ifferentiating Equation (14) 



(29) 



(30) 



^SV A >>' 

^ - r^?7^ 

Similarly, differentiating Equation (16) 






AT, A/, 






# 

w> 






> ( 31 ) 
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But now the values of and cannot be evaluated to 



the first order from the linearized expressions for and 
in Equations (13) and (17) respectively because the q^ term has 
vanished in the process of linearizing these equations. Therefore 
the differentiation must be carried out upon the exact expressions, 
Equations (10) and (8), and the results then linearized. 

From Equation (10) 




where 






and the time derivative is 




Therefore 











£09 & 









( 33 ) 
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and upon differentiation 



which, by retaining only first order terms may be written 

A£y ^ ^ y.ACfy-f,l 

Similarly 

Xy 

The various derivatives of £> must be computed from 
Equation (12) 




Thus 



/ 





caCfy-fs) ^..rYfy^s') 



and approximately 










( 35 ) 



Now 



y C 



C40J <S> 
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and substituting the results of Equation (35) and linearizing 

Similarly, from the results of Equations (35) and (13) 

fs'' ^ ^ 

Substituting Equations (34), (35), (36), and (37) into Equa- 
tion (33) and retaining only first order terms 






Jv 






or 



To find the same procedure is followed. Thus, from 

Equation (8) 



Ts 









and the time derivative is 












• I 
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Therefore 

^^3 _ /V/ a//7 cP <9 ^ « 

^f3 

. ^ <?«»;» <9 

ts'/j 

/^/ /^ Jc.s<S> /Y^-c <S> /^ ^ <g (40) 

Substituting Equations (34), (35), (36), and (37) into Equation (39) 
and retaining only first order terms 



^ ^ ^ Cf,f. /J + (fY -fj) 

_ . ^ /yj 



-^ - ^^f>CJ,-A) 7sU-A) /4JX 

The second term of the Lagrangian can now be written by 
substituting Equations (30), (31), (32), (38), (39) and (41) into Equa- 
tion (30) 

4f, = 



-26- 



- ''fJ 

^ '^//* ' '^/' - rMr ^'^J 



or 



iP7 ~ _ /%^j. 









/%/^l 









(42) 



The generalized force is 



= -W or^ j 

V- 



Therefore, referring to Figure 5, 






Substituting Equations (3), (9) • and (10) to eliminate the extraneous 



coordinates . and q^ , and collecting terms 




Figure 5. Schematic diagram of external forces acting on system. 



Carrying out the indicated differentiation with respect to q 
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Subatituting Equations (34), (36), and (37) into this expression, 
and using the approximations fy ^ J-* ^ ^ 

becomes 



Finally, retaining only first order terms, the expression for the 
generalized force reduces to 

Therefore, by substituting the results of Equations (28), 
(42), and (43) into Equation (18), the equation of motion in the q, 
direction becomes 






(V, 
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DCRIVATION CF EQUATION OF MOTION 
FOR COORDINATE q4 

The first term of the L>agrangian expression is obtained in 
an identical manner to that employed for coordinates * 

and q^ . 

From Equation (19) 




and, differentiating Equation (14) 




so that 




Similarly, differentiating Equation (16) 




so that 






Therefore , 




or 
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and the first term of the Lagrangian expression becomes 



From Equation (19) 



A. 



i - -1^) 

Lfferentiating Equation (14) 



so that 

Ajkr - « 

Similarly, differentiating Equation (16) 

f -P- 

so that 

^ 

^ 

Here again the values of - ^ =- and cannot be 

ated to the first order from the linearised expressions for q 



(45) 



(46) 



> ( 47 ) 



^ (48) 



evalu- 

and 

(m 



q^ in Equation (15) and (17) respectively because the q^ term has 
vanished in the process of linearizing these equations. Therefore 



the differentiation must be carried out upon the exact expressions. 
Equations (10) and (8). and the results then linearized. 



The time derivative of , defined by Equation (10), has 
already been found from Equation (32) 












Therefore 



^ fjyj, <s>l J> 







<0 y Afz 






,fjy 




A c , 




The value of 


follows 


from Equation (ll) 






— ees Cpr -J'j) 


and upon differentiation 















(49) 



which, by retaining only first order terms, may be written 
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Similarly 

^ (^fy -fs) 

Jr 

The various derivatives of d> must be computed from 
Equation (12) 







Thus 



/ ^ 






J. 






•iff 









<?& _ _ Jx. JrroX'f^-/^) -JJ 

^fr " JJ-Jj^'-'TYrr-p) 

and approximately 



^ ^ _ J 

^ JJ Jr 

Now 

_ (Sbfj<S> J<S> 



( 51 ) 



and, substituting the results of Equation (51) and linearizing, 









J 



Similarly, from the results of Equations (51) and (13), 



(52) 



^ s= _ <S ^<£> ■= J ^ J ~-J 

7^ {a-' 



( 53 ) 




I 
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Substituting Equations (50), (51), (52), and (53) into Equa- 
tion (49) and retaining only first order terms 



To find the same procedure is followed. The time 

derivative of q^ , defined by Equation (8), has already been found 
from Equation (39): 

Therefore 



jiA= j/a (-^y^<o)e , A^> S ^ 



. <?» J <5 _ ^ ^ 



<ro^<9 









X <g 



(55) 
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Substituting Equationa (50), (51), (52), and (53) into Equation (55) 
and retaining only first order terms 



The second term of the Lagrangian can now be written by 
substituting Equations (32), (39) « (47), (54), and (56) into Equa- 
tion (46) 






^r _ _ 

^ ^ f' ‘ V* ■*/^/ 






A/, -f /V* 0 ^ 



( 57 ) 
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The generalized force is 

^ Work^ 

Therefore, referring to Figure 5, 



^ -A ^ -4 ^ p] 

- 

Substituting Equations (8), (9) • ^nd (10) into this expression to 
eliminate the extraneous coordinates . and q^ 

■>- ff - ^ «> <£> - U*Ji) --"/SZ 

and collecting terms and performing the differentiation 

Substituting Equations (13), (34), (36), and (37) into this expression, 
and using the approximations and becomes 



-36- 



* f''- 

Finally, retaining only first order terms, the expression for the 
generalized force reduces to 

Therefore, by substituting the results of Equations (45), 
(5 7), and (58) into Equation (18), the equation of motion In the q^ 
direction becomes 



('-p - ^ - 






( 59 ) 



Thus Equations (23), (27), (44) and (59) are the four La- 
grangian equations of motion corresponding to the four principal 
coordinates of the system. 
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ANALYSIS OF THE EXTERNAL FORCES AND REDUCTION 
OF THE EQUATIONS TO DIMENSIONLESS FORM 

In carrying out the stability analysis the state of vertical 
motion corresponding to the direction will be considered con- 
stant during the perturbation. That is, the values of and q^ 
will be fixed, and their magnitudes will be determined by the 
time -instant at which the stability of the system is investigated. 
Therefore the equation of motion in the q^ direction, Equation 
(27). will not be considered further In this analysis. 

The thrust force, F, of the booster may be expressed in 
terms of the mass rate of flow, M^ , and the specific impulse of 
the fuel, Ijf, , as 
/" = - 

The negative sign is necessary because M^ is numerically a 
negative quantity as it occurs in the Lagrangian derivation. Thus 
the thrust, F, will become a positive quantity when evaluated nu- 
merically in the next section. 

Also, it was found in the first p>art of the derivation that 
the jet damping force, D, is proportional to the quantity 
In order to make this a numerically positive quantity, it will be 
written as 

2 ? - - 



38- 



The moment o£ inertia, , can be written as the product 
of the mass, , and the square of the radius of gyration, k^. , 
that is 



The time derivative of the moment of inertia, , is 



J , Howe ver , the second term is 



small compared to the first term and will be neglected in this 
analysis . Therefore the following approximation can be made 

Making these four substitutions in the equations of motion 
in the q^ , q , and q^ directions the following equations are ob- 
tained from Equations (23), (44), and (59) » respectively. 

From Equation (23) 



From Equation (44) 









/■- 
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{/r.*/^^}-^ * ATf^^^x ^ 'V/^ 

From Equation (59) 



( 61 ) 



/vj/% /!> >i /*)?*/-5>i ^ A ^ 









' AZ-f^/^i. 






(aZ *Afx)^ 



4A 



o 



(62) 



To write these last three equations in dimensionless form 
introduce the dimensionless coordinates . Q, , and Q^ , and the 



characteristic units of time = I,^ , velocity s 


lj;.g* acceleration 


s g* 










Thus 










‘JJ = -^ 


=-A- < 

^s,<y 








II 


II 






> (63) 






II 
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Making the above substitutions into Equation (60) and divid- 
ing the resulting equation by M, g, the following dimensionless equa- 
tion is obtained 









■J 1:7^ 



(0 f 






A 



4 



A±. ^ 

r.;j 



/* 75 , Tr^ 
A/fy/Vi 






/ /A Ar ^ 

z;<j /^z 



/• 



/^zTz^ 

/T’* 









(64) 



Similarly, substituting the values of Equation (63) into 
Equation (61) and multiplying the resulting equation by 7^'* 
the second dimensionless equation is obtained 
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Finally, substituting the values of Equation (63) into Equa- 
tion (62) and multiplying the resulting equation by the 

third dimensionless equation is obtained 



Affl' /ft ^ * 

_Am. ^ A 3 fAtAf] -A. A A 
M* A _ Ar ff _/%, ) 



( 66 ) 



Equations (64), (65), and (66) are the three dimensionless 
equations of motion of the system and contain the following eight 
dimensionless parameters 





A 


/ft. 


A 




^y9 




A 


A 


At 




At 








AT 
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III. NUMERICAL INVESTIGATION OF STABILITY 



Rather than attempt to find a solution for Equations (64), 
(65), and (66), it was decided to introduce selected values of the 
first six dimensionless parameters and, treating and 

as variables, to use Routh's Stability Criteria to determine 
the range of combinations of these two variables in which the 
system would be stable . 

To simplify the writing of the determinantal eq[uation 
necessary for the use of the Routh criteria, introduce the fol- 
lowing notation where 



A^^. coefficient of in Equation j 

coefficient of in Equation j 

C^y coefficient of Qy in Equation j 

and 

j = 1 corresponds to Equation (64) 

j = 2 corresponds to Equation (65) 

j = 3 corresponds to Equation (66) 

With this notation Equations (64), (65), and (66) may now 



be written as 
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C„ C,+ B.. Q. + C3 + A3. O, + + A,. = 0 

B,i C. -f Q, + B,, Q, + A,^ + B^^ + A^^ = 0 

Bj, 0.+ C^3 Q, + B3, C, + A,, O, + C„ Qv + B^, d\ + A^, Q., = 0 



where, letting ^ and^= ^ , the coefficients in Equation 

(67) are 



A3, 


3 




A3. 


= 




Ai. 




/y. ( 


Ay, 


= 


A/i. Xs^ f ^ Xry> )^ ^ 

/vr^/^. r 7:^ ^ 


A.. 


S 




Ay, 


3 




B« 


S 






S 


- A, A* ) 

Af,-t-ATz 1 /^/ / 


B.s 


S 


/V: 

/<fi 1 f^- J 
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Af.*Aft. ^o/tj ' 



/V, T,J^<^ ^ 



A , /r^\ ^r- ^ 

_ / ^ /^ I -^j- , /V< -^/» yif7 

/ ^.*^A Jzp^ /j 

aTkI z^*^ 

" T^zV-^T^J 7^ 

h^) 

^lti‘(:-S);i^^jh-^ 



4z 

5^ 






-_f^l // 



_^ //^ ^y_L 
j^/^r /^/^ 



^y_L^ 



'/^ 

mj-o 
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Inasmuch as oscillatory solutions are of main interest. 



assume 






Then 


4 - 




and 







For the system of homogeneous differential equations (67) 
to have other than a trivial solution, the determinant of coefficients 
of must vanish. That is 



C„ A + B„ A A + A,, B^^A + 

B,,. A Of X A + Bj^ A + Qj,* A + B^t A + A^ 

B,j A Cjj \ + Bgt A + Ajj C^s ^ + By, A + 



■■ 0 (67a) 



Expanding the determinant of (67a) 

C„ FA"^+ (C„ H + B„ F)A%- (B„ H + C„ J - B.^S + B„\V) x"' 
t (C„ K + B„ J - B,^T + B,,X) A^ + (C„ R + B„ K - B,;^U + B.,Y) A* 
+ (B„ R - B.jV + B„2) A = 0 
which can be written 

P., a"' + a" + A^ + 15^ A + =0 (68) 

where 

P, = C.. F 

p_ a C„ H + B„ F 

p^ = B„ H + C„ J - B.^S + B„W 

p^ = C„ K + B„ J - B,j T + B„ X 



J 

I 




i 

1 

b 




T 
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= C„ R + B„ K - U + B„ Y 
= E„ R - B,*V + B„Z 

and 



F 

H 

J 

K 

R 

S 

T 

U 

V 

w 

X 

Y 

z 









“ss'-'yi 

Vj ” ^KZ^JS ~ ~ ^S3 



= B,, - B,, C„ 

“ ■*■ B^^ Byj ~ A^^C 



S3 *~^33 "V/- 



“ ^3/ ^ys ~ A»j A^^ 

” Bj,C^ - Cj;j 

“ A,. C + E, B , - A C , 



'SZ ^3Z 



“ Aj/ A^^ - Ay^ A ^ 



The system will be stable If the coefficients of Equation (68) 
and the following terms of Routh's stability criteria are all of the 
same sign (Reference 3) 






r. 



/^3 - 



> ( 69 ) 



- 

^ -/? y^s-U 



3 . 
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NUMERIC AL EXAMPLE FOR ' V-2" TYPE ROCKET 



To find a numerical solution to the problem the "V-2" rocket 
(weight approximately 28,000 lbs.) was chosen as the main rocket, 
Mj, and a small five-second booster providing an initial accelera- 
tion of 1-1/2 g was selected as the booster, M2. 

The dimensionless parameters of the system were then es- 
timated to be 





-p- = 




4-V 


4- ^ 7.7 s- 








= <7 (at launching) 



The stability of the system was then to be investigated for various 
values of the two variables and . 

Because of the inherent complexity of the constants of Equa- 
tion (69), a vast amount of time would be required to investigate 
every likely combination of the variables and . How- 

ever, the writer optimistically set out to evaluate these constants 
for various combinations of a and^<^ . Inasmuch as no combination 
was found where even the constants p^ to p^ were of the same sign, 
much less where the test functions of Equation (69) were of the same 
sign, it was decided to reduce the complexity of the system in order 



to expedite the computation. 



I 
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Accordingly, and Q, were assumed to be zero, that is 



the system was restricted to vertical and rotational movement of 
the center of gravity. This assumption reduced the system of equa- 
tions to two, namely, Equations (65) and (66), and reduced the labor 
Involved in applying the Routh stability criteria substantially. Once 
a stable range of a and^ had been determined, these values were 
then to have been used in the general case for closer investigation. 

STABILITY OF SYSTEM WITH TWO DEGREES OF FREEDOM 

If both Q, and Q, are assumed to be zero. Equations (65) 
and (66) may be re-written (using the abbreviated notation intro- 
duced earlier) as 




and assuming as in the more general case, the determi- 



nant for this case becomes 







A + A + Aj 



The expanded determinant may then be written 
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■*■ ■♦■ AjjjC^y + Bjj^ “ ^K-z 

'*' ■*■ ~ ®<rz^jj *■ A^^Bj^)A + - A^^Aj^) = 0 (70) 



Furthermore, since the parameter — ^ — 
the last term of Equation (70) becomes 



equals zero at launching, 




Therefore Equation (70) reduces to a third order equation of the 
form 

/? -f- -h * /^3 ^ 



where 



P. “ ^i3t.^irs ~ 

I> = “*■ ~ 

P( ~ ^3A^vs ®a* “ ^Z ^JS ~ \z^fr ~ 

Pj ' '*' ^Jl^YS “ \z^3 3 

and the Routh stability criteria require that • Pj » 

f/l - -^^3^ all be of the same sign. (Reference 3). 

The coefficients of Equation (71) may now be evaluated by 



substituting the parameters previously selected for the "V-2" rocket 
with five second booster. Collecting terms these coefficients may 
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be expressed as follows 
where 



=• ^zV * /«3' 



^ = Z. 7 :T.^ o' 



= <-J 



(T2) 



(73) 



where 

^ = -/.<^/ . 

d, = S,ZLz-i^^ d./d-nyo '"^ Y.Z . ya'^ 

= Ctg tC^ * c 6 ■*■ C:x. (74) 

where 

4 /iP~"y^^ 

^ = J.d-3 »/C>~'^^i- 'y.73 • AZ/ ‘ /a'" 

d. = ~f/.7<i • /o'^*^^,77»/o'^ * £ J7./q 't. fT 3 >/o'^ 

~ Sjo' -f- ^3 (75) 

where 

=: 4 yo~'^^^ 

= -^d.a<^./d>''''/f*-^ S^/Z ’/o'^) 

= 7.o3>/o'/f’* S./Z‘/o'^ * /./^»/0~^ 7./^‘/o 
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Before computing numerical values of the coefficients of 
Equation (71) it will be profitable to draw certain conclusions from 
an inspection of Equations (72) through (75). 

a is defined as the quantity ^ , and ^ is defined as the 
quantity , where ^ and are inherently positive quan- 

tities by the physical nature of the problem. Therefore, it can 
be concluded that for any real solution of the physical problem, 
both a and^ must be of the same sign, that is, either both posi- 
tive, or both negative* 

In Equation (72) the sign of p^ is independent of and 
takes the sign of a. Therefore numerical values of p^ need not 
be determined at this stage of the investigation since the sign of 
p^ will be evident by inspection. 

In Equation (73) a, is always negative, and c, has roots 
at ^ = 0.456 and = -1*93, 



Since both a, and c, are 
negative for /<f < 1*93 and 
since a is negative when 
is negative , it follows 
that P| is always posi- 
tive in this region. 

Negative a yields negative p^ 




from Equation (72). Therefore p^ 



and p^ will be of opposite sign for all values oi /S < 1.93, and. 
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by Routh's stability criteria, stability is impossible for this case. 
By analogous reasoning stability is impossible in the range 
0 <^< 0.456. 




In Equation (74) a^ = 0 when = 0, and a^ is negative for 
all other values ol , Furthermore, b» has a root at =2.54x10 
and c* has a root at 
- -2.68x10^ . Since 
a^ , bf , and are all 
negative in the range 

it can be concluded that stability is impossible in this range by 
reasoning analogous to that in the case of p^ . 

In Equation (75) a^, = 0 when = 0, and a j is positive for 
all other values of /3 . 

Furthermore b, has a 
root at/? = -1.34x10 
and at = 0 . Also 
c 3 has roots at /^ = 

0.456, = -1.92, 

and /<? = -2.68x10'^. 




No conclusions can be 
drawn in this case . 

To summarize the conclusions of the preceding paragraphs. 
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the possibility of stability need only be investigated in the ranges 
- 1.93 </^< 0 and 2.54xl0''^<^ < 0.456. 

In evaluating the coefficients » P* • P? follow- 
ing procedure was adopted. Selected values of within the 
positive range 2.54x10 < 0.456 were used to find the con- 

stant coefficients of Equations (73), (74), and (75). Then the 
positive roots of the resulting equations were computed as func- 
tions of a. Negative roots were not investigated since a must 
be of the same sign as . These values are tabulated in Table 
I. 

Similarly, values of the coefficients and corresponding 
negative roots of Equations (73), (74), and (75) were computed 
and are tabulated in Table II. 

The data of Tables I and Table II are plotted qualitatively 
in Figure 6 and Figure 7. 

In Figure 8 and Figure 9 are plotted the logarithms of 
maximum and minimum limits of a for which p, , p^ . and p^ 
are each positive. 

It is apparent from Figure 8 that there is no single value 
of Q for which p, , p^ , and p^ can all be positive since 
for positive p^ greatly exceeds for positive p^ for all 

values of 2.54x10 </ 3 < 0.456. 

Likewise, in Figure 9 there is no single value of a for 

m'. 



I 
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which • Pt • P 3 can all be negative since la^,,l for nega- 
tive greatly exceeds for negative p for all values of 

-1.93 </^ < 0. 

Therefore, it is clear that .'Equations (65) and ( 66 ) do not 
have a stable solution for any combination of a and<<^' for the 
parameters chosen. 

STABILITY CF SYSTEM WITH REDUCED MASS RATIO 

In view of the negative results obtained in the previous 
case it was decided to investigate the effect of changing the par- 
ameter . The following set of dimensionless parameters 
was then selected 



_ /.s- 7^*1^ _ -z.i/ 



A 

A 



/. S' 






f*-. = <D 



(at launching) 



Note that was decreased from 7.75 to 1.5, ^ 

was increased from 1.96x10 to 6.72x10 and was In- 

9 

creased from 2.48x10 to 8.75x10*. The remaining parameters 
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were unchanged. Thus, this second set of parameters describe 
a system essentially the same as that first investigated except 
that the mass ratio has been markedly reduced. 

/Vi 

The coefficients of Equation (71) are again evaluated by 
substituting these new parameters. Collecting terms the coef- 
ficients may be expressed as follows 



p^ = ^ ^ 



(76) 



where 5 : <?<•/<? 










(77) 



where « /o 

_ 



<3^ c/ + Cz 



(78) 



where a^ = . z.zz ^ /o 

^ . zo z. c. /o' 

Q S - fz.Xh,/o''^\ ^ V./9*/0 ^^ 

^ ‘ /3 



^ V./Z^/o 



/<■/» /<? 



-T 



CTj o(. ■+ Cj 



( 79 ) 



where a_, * 7. 3z t /o~'“/<3*' 







V 



- 7. sz* 7o''/3^*S.ZZ “70'^ - . /p 
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As previously stated, a and must be of the same sign. 
Again, the sign of is independent of y and takes the sign of a. 

In Equation (77) a, 
is always positive. Also 
c, has roots at /S =0.68 
and at /? = -2.84. 

In Equation (78) 
a, = 0 when ^ ^ c , and 
is negative for all 
other values of . 

Furthermore b* has roots 
at /^= -5.6x10^ and at 
= 2.55x10 Also 
has a root at/^= -1.12x10'^. 

Since a^ , b^ , and c, are all negative in the range < 2.55x10 

stability is not possible in this range. 

In Equation (79) a^ = 0 when /f = 0, and is positive for 
all other values of . Furth- 
ermore b^ has roots at/^ = 0 
and at /^ = -5 .62x10^ . Also, 
has roots at /<^ = 0.625, 






/<? = -2.83 and/5 = -1.12x10 . 
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It may be concluded that the possibility of stability need 



As in the previous case selected values of both positive 
and negative/^ were substituted into Equations (77), (78), and 
(79), to compute the constant coefficients of Equation (71), and 
the resulting equations were then solved for maximum and min- 
imum allowable values of a. The results are tabulated in Table 
III and Table IV, and are plotted qualitatively in Figure 10 and 
Figure ll. 

Figure 12 and Figure 13 are again logarithmic plots of 
the maximum and minimum limits of a versus the logarithm 
of ^ . 

In Figure 12 it can be seen that stability is possible in 
the range 10 < < 4.5x10^ only if for positive p^ is 

greater than <x„,^ for positive p^ . (The curves appear to co- 
incide.) However, an analysis of Table III reveals that for all 
values of 10*'^^< 0,68 both c, and c, are negative whereas 
is negative while a^ is positive in this range. Therefore 



Thus it follows that is greater than in this range 

and stability is impossible. 



only be investigated outside the range 2.55x10 , 




while 




Similarly, in Figure 13 stability is only possible if 
for negative p^ is greater than \a for negative p^ , 







9 



1 
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However, an analysis of Table IV reveals that for all values of 
-2.8 -10’ both and Cy are positive while a^ is negative 
and aj is positive in this range. Therefore p y y 

while and thus is greater 

tl^an . 

Therefore it is only necessary to investigate more close- 
ly the range 0,68 <'/^< 4,5x10^ and -2.8. 

Returning to Equation (71) and multiplying the analytic 
expressions for the coefficients of p^ and p^ the following equa- 
tions are derived if only first order terms are retained. 









/^ = Cty a(^ -i- Sy 0L Cy 



( 80 ) 
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Factoring ^ from iCquatlon (81) 









I /y./zYf ^ ^J ('~ A*' ^ 










(82) 



Comparing Equation (80) with Equation (82) it can be seen 
that both a^ and b_, are multiples of a^ and b^ respectively by the 
factor <^1 . For positive/^ , stability exists only if 

is greater than . that is , if ^ ^ 




m 



-60- 



But the factor 

/*)f, ^ A-/t_ 

ically negative ao that the condition for stability can be 





is numer- 
written 



Substituting the values of c^ and c, from Equation (80) and Equa- 
tion (82) this condition becomes 









which, upon collecting terms, becomes 






-?*/• ' 



Inserting the numerical values of the parameters, the second term 
of the* denominator is of the order 10 since - ~ C.7i ' 

Therefore this term can be neglected in comparison with the first 
term of the denominator and the expression simplifies to 









^S- 



( 83 ) 



Solving this inequality with the parameters chosen, the result is 



^ 1.62 



( 




A.a an example, chooce ^ s 2. Inserting values from 
Fable III into Equation (78) and Equation (79). for = 2 the 
following values of a are obtained 




" /o'^ ^/ - Z. /0~^) 

c^3 5r - /a^f / ~ 

Thus changes from positive to negative when 
Oi = 5.6x10^(1-2.54x10^), while changes from negative to 
positive when Cj = 3.6x10*^(1-9.83x10''^, and both ]s and arc 
positive lathe range (5.6x10 -5 .52)< o< (5 .6x10^-1.42) . Re- 
sults with a similar sensitivity are obtained for larger values 
of /3 but smaller values of a. (It can be seen from Figure 12 
that the product cf a and ^ is approximately constant and equal 
to lO"^.) 

There is yet to bo applied the final Routh criterion, 
namely that /? - be positive also. Choosing at its 

smallest allowable value, say 0^, and calling the correspond- 
ing value of some small positive value, say r, ^ . 

Thus all the stability criteria are satisfied and the system is 
stable. Therefore, it can be concluded that the system is stable 



-62- 



for > 1.62 but that, for any given value oi/f in this stable 
range, a is extremely restricted. 

Finally, there remains the negative range of^ to be 
investigated, namely -10^</^< -2,8 . But, for negative values 
ol ^ must be greater than -1.62 by Equation (83) in order 

that and both be negative. Therefore stability is not 

possible for negative values of 



0 ( 
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IV. DISCUSSION CF RESULTS 

Only two possible sets of parameters for a launching sys- 
tem of this type have been investigated. In the first case where 
the mass ratio was 7.75 it was found that stability was impossible, 
whereas, in the second case where the mass ratio was 1.5 it was 
found that stability was theoretically possible for values of be- 
tween 1.62 and 4.5x10^ if the effect of lateral motion of the center 
of gravity of the system was neglected. For any given value of ^ 
in this range stability was possible for only one corresponding 
value of a . 

Returning to the physical aspects of the problem consider 
now the meaning of the variables a and ^ ^ is defined as 

where is the distance from the center of gravity of the booster, 
M* , to the pin connecting It to the strut to the main rocket, M, , 
and kj is the radius of gyration of the booster, . For any sys- 
tem of this type k^ will always be of the order of one foot or 
greater. Therefore, will always be equal to or greater 

than , That is, as a first approximation, ^ is approximately 
numerically equal to . Thus, from a practical consideration 

, and therefore ^ , must certainly not exceed some small num- 
ber of the order of 1 to 10, depending on the size of the booster, 
for otherwise the pin will be located too near the booster exhaust. 
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Now a is defined where X is the length of the strut 

from the pin in the booster to the center of gravity of the 
main rocket^ M, . Here again there is a practical upper limit 
for X, , first, because the strut was assumed to be of zero mass 
in the derivation of the equations of motion, and second, because 
the strut was assumed to be of infinite rigidity. A long strut 
would present structural problems in rigidity since it is essen- 
tially a column acting under both axial and bending loads. Re- 
stricting Xt. to values of the order of a few feet it is clear that a 
should not exceed 100 or 200 at most, since these values corres- 
pond to an Xt of 25 to 100 feet. 

In view of the foregoing remarks consider the product of 
the two variables a and . Thus . If 

is of the order of 1 to 10 feet, say, and X, is restricted to 25 to 
100 feet, then it can be seen that the product should not ex- 
ceed 25. Returning to Figure 12 note that, in the range 1.62 <^< 
4.54x10^ , a is approximately constant and equal to about 10 
Thus, even though theoretically stable, this system is of no prac- 
tical significance for the mass ratio investigated. 
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V. CONCLUSIONS 

1 , A compound pendulum system of launching stabilization is not 
stable for a mass ratio 7,75, 

2, Such a system is theoretically stable for a mass ratio 1.5 if 
the effect of lateral motion is neglected. However, even in this 
restricted case the system is of no practical significance due to 
the excessive length requirements for the strut from the main 
rocket to the booster, 

3, Since the system is quite sensitive to the particular param- 
eters chosen, several more investigations of specific cases are 
required before any general conclusions can be drawn as to the 
practicality of this method of stabilization. 



Data from ii^quations (73), (74), and (75) 
(2.54xl0'<^< .455) 
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Data from Equations (73), (74), and (75) 

(-1.91 < 0) 
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TADI-E III 

Data from Equations (77), (78), and (79) 
(^>2.55 X lu"^) 
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Figure 6, Variation of p » and p with a for several values 
of positive ^ , ( ^ “ 7.7S) 





Figure 7. Variation of * P* » P 5 with a for several values 
of negative . ( -^ = 7,75) 
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Figure 8. Minimum and maximum limits of a which yield posi- 
tive p for positive^ . ( ^ = 7o75) 
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Figure 9. Minimum and maximum limits of a which yield nega- 
tive for negative 7.75) 



-74- 









Figure 10, Variation of p ■ * and with a for several values 

of positive ^ ^ - 1,5) 
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Figure 12o Minimum and maximum limits of a which yield positive 
p. for positive ^ ^ sr 1,5} 
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Figure 13 o Minimum and maximum limits of a which yield nega- 
tive i> for negative ^ ^ 
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Figure 14. Schematic diagram of system showing inertia 

forces for use in Newtonian derivation of equa- 
tions of motion. 
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APPfiNDIX 



NEWTONIAN DERIVATION OF EQUATIONS OF MOTION 

The equations of motion for the two dimensional system 
may be derived by applying Newton's Second Law of Motion pro- 
vided that proper consideration is given to the inertia forces of 
the system. In Figure 14 these forces are indicated by broken 
lines. In this derivation the main rocket, M; • and the booster, 
Mi , will be treated separately as free bodies and the reaction 
of the pin connection will be replaced in Figure 14 by its two 
components, f, in the vertical direction, and f* in the horizon- 
tal direction. The pin is assumed to be frictionless and there- 
fore no moment is transmitted. 

Treating the main rocket, M , , as a free body first, the 
sum of the vertical forces acting on it are 



Similarly, the sum of the horizontal forces acting on M, are 






and therefore 




( 1 ) 



zH =4 



and therefore 



A ^ 



( 2 ) 



-80- 



Finally, the sum of the moments acting on M, are 
or 

( 3 ) 

Treating the booster, , as a free body, the sum of the 
vertical forces acting on it are 
TV ^ -P^yyf 

and therefore 

Furthermore, the sum of the horizontal forces acting on are 
and therefore 

4 . - ^ ^ (5) 

Finally, the sum of the moments acting on are 

TM=i^f^ t i^p V- ^ pfA-^A) 

or 

( 6 ) 

Substituting the results of Equation (1) into Equation (4) to 
eliminate f, 

Similarly, substituting the results of Equation (2) into Equation (5) 
to eliminate f^ 



r /V * 









( 8 ) 
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Substituting the results of liquations (4) and (5) into Equation (3) 
to eliminate £, and 






(9) 



Similarly, substituting the results of Equations (4) and (5) into 



Equation (6) to eliminate f, and f^ 









( 10 ) 



Assuming the angles q^ and sufficiently small such that 
terms of order q* may be neglected, the following approximations 
can be made 

Therefore, Equations (7), (8), (9)» and (10) may be rewritten in 
simpler form. Thus 
From Equation (7) 



f^- 






From Equation (8) 
From Equation (9) 






(13) 



] 




I 




I 



I 
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From Elquation (10) 

p -^pj,]4jy (14) 

The following approximate relations between the principal 
coordinates , q^ , q^ , and q^ , and the extraneous coordinates 





* / 




(15) 








(16) 


% 


* f 




(17) 








(18) 



Taking the first and second time derivatives of Equation (15) 



cing the first and second tlm< 



Taking the first and second time derivatives of Equation (16) 
/'>’ _ ) ( 21 ) 



/; -/i - <«) 

Taking the first and second time derivatives of Equation (17) 
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( 23 ) 



and 





(24) 



Taking the first and second time derivatives of ^^Iquation (18) 



liquation (11) can now be rewritten in terms of the principal 
coordinates only by substituting the results of iSquations (22) » (25) 
and (26) for the extraneous coordinates 



Note that Equation (27) is identical with Equation (27) of Part II of 
this paper. 

In a similar manner, Equation (12) can be rewritten in terms 
of the principal coordinates by substituting the results of Equations 
(20), (23), and (24) for the extraneous coordinates 



which is identical with Equation (23) of Part II of this paper. 

Substitution of Equations (23), (24), (25), and (26) in Equa- 




(25) 



and 




(26) 





tlon (13) yields 







1 
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= - 7^ f ■^^p-p * -^J '”> 



which is idantical with Equation (44) of Part II of this paper. 

Finally, substitution of Equations (23)* (24)* (25)* and (26) 
in Equation (14) yields 









( 30 ) 



which is identical with Equation (59) of Part II of this paper. 
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TABLiC OF SYMBOLS 




D - Jet damping force 
F - Thrust force 
g - Acceleration of gravity 

I - Moment of inertia 

- Specific Impulse 
k - Radius of gyration 
M - Mass 

M - Mass rate of flow 

q > Lagrangian generalized space coordinate 
Q - Dimensionless Lagrangian generalized space coordinate 
^ - Lagrangian generalized force 

T - Kinetic energy 



- 86 - 



PviFiSRfCNC JS 



1. J. M. J. Kooy and J. V/. H. Uytenbogaart; "Ballistics of 

the Future", (1946), Chapter XI. 

2. Theodore von Karman and Maurice A. Biot: "Mathematical 

Methods in Engineering", (1940), Chapters III and VI. 

3. E. J. Routh: "Advanced Rigid Dynamics", (1905), Vol. II, 

Article 302. 



JUL 2 



B t NOCtl'V i 



15648 



Thesis 

N9 

An analysis 
xil-om rocket 



Horris 

of a conpoiiiKi pend- 

suspension. 



lS 



jn 



Library 

S. Naval Postgraduate Sckool 
Monterey, Caliloniia 



